Spectroscopic Evidence for the Localization of Skyrmions near zv=l as T — i> 
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Optically pumped nuclear magnetic resonance measurements of ^^ Ga spectra were carried out in an 
n-doped GaAs/Alo.iGao.gAs multiple quantum well sample near the integer quantum Hall ground 
state v=l. As the temperature is lowered (down to T fs 0.3 K), a "tilted plateau" emerges in 
the Knight shift data, which is a novel experimental signature of quasiparticle localization. The 
dependence of the spectra on both T and v suggests that the localization is a collective process. 
The frozen limit spectra appear to rule out a 2D lattice of conventional Skyrmions. 



One of the most surprising twists in the recent history 
of the quantum Hall effects [|l| was the prediction ||] that 
novel spin textures called Skyrmions can be the charged 
quasiparticles introduced by small deviations (| 5 z^|) from 
ferromagnetic quantum Hall ground states g] (e.g., at 
Landau level filling factor v=\ or i). A Skyrmion has 
an effective number of spin reversals K and "size" A that 
are determined by the competition between the Coulomb 
energy (which increases both) and the Zeeman energy 
(which reduces both). Qualitatively, this cylindrically 
symmetric spin texture has a down spin at r=0 and a 
smooth radial transition to up spins at r=oo. In be- 
tween, the nonzero XY spin components have a vortical 
configuration [Hj]. The addition of Skyrmions to the 
v=\ ground state was predicted to result in a rapid drop 
in the electron spin polarization, as \^v\ is increased [p|. 
Several experiments are consistent with this P-pj and 
other predictions |pHl3[ of the Skyrmion model. 

These developments stimulated many theoretical stud- 
ies of Skyrmions. One of the central questions that 
emerged was the nature of the many-Skyrmion ground 
state. Would Skyrmions form a crystal [p^-[20|, and, if 
so, what symmetries would it possess? Does disorder |2|] 
affect the Skyrmion size A and spin number ii" as T ^ 0? 
While these questions are still under active investigation, 
basic aspects of the Skyrmion model have not yet been 
tested experimentally. For example, the detailed shape of 
a Skyrmion has not yet been measured with local probes, 
presumably because the quasiparticles are delocalized at 
high T. 

In this Letter, we report the first spectroscopic evi- 
dence for Skyrmion localization. The multiple quantum 
well sample used in this work was previously studied at 
higher temperatures ^,|l^ . The new data presented here 
were obtained by extending the optically pumped nuclear 
magnetic resonance (OPNMR) technique |23] to lower 
temperatures (T « 0.3 K) as described elsewhere p3| , p4[ . 

Figure H shows some OPNMR spectra for v close to 
one. Nuclei within the quantum wells are coupled to 
the spins of the two-dimensional electron system via the 
isotropic Fermi contact interaction Ea], which shifts the 
corresponding well resonance (labeled 'W" on Fig. ITIc) 
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FIG. 1. Several ^^Ga OPNMR emission spectra (a-c, solid 
lines, offset for clarity), acquired at low T for different v, 
in Btot=7.03 T (/o=91.36 MHz). The same 'W" resonance 
is obtained using either conventional NMR absorption (a, 
dashed line) or OPNMR emission (a, solid line). The dot- 
ted lines (a-c) are fits to W and "_B" , which assume identical 
quantum wells (see text) . 



relative to the signal from the barriers ("i3") RM. We 
define the Knight shift Kg as the peak-to-peak splitting 
between W and B. The third resonance (labeled "AT" 
in Fig. ^) may be due to nuclei near "defect" sites. By 
scaling equilibrium signal intensities, we estimate that X 
is due to only '--^ 2% of quantum well nuclei ||26| ; for now 
we ignore it. 

For v—\^ all spectra (e.g.. Fig. y:) are well described 
by the same two-parameter fit (dotted lines) ||2^,|2J] that 
was previously used for all spectra at v—\^. Note that 
this fitting function has no explicit dependence on the 
(x,y) position of nuclei along the quantum well, despite 
the fact that NMR is a local probe. This is because the 
fit is generated under the assumption that all electron 
spins are delocalized, so that {Sziy, T)), averaged over the 
NMR time scale (~20/isec), appears spatially homoge- 
neous along the plane of the wells, and thus the resulting 
lineshape is "motionally- narrowed" pq |. In this limit, 
measurement of A'5 reveals the "global" , time-averaged 
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FIG. 2. Ks(u) (solid symbols) and Tt^{u) (small open sym- 
bols) at (a) r=3.7 K, (b) T=l.& K, and (c) r=0.4 K. The 
large open symbols in (b) are Ks data previously reported 
for this sample. Here, r„ is the full width at half maxi- 
mum (FWHM) for the "W" resonance. The filling factor 
{v—nhc/ (eBtotcosO)) is varied in situ by tilting the sample 
(0°<6I<37°). In practice, all the low temperature (T<1.6 K) 
data were acquired following a cooldown from at least 1.6 K 
at fixed V. The 2D electron density (n = 1.52 x lO^'^ cm~'^) is 
inferred from the Ksiy) peak in (a). Solid lines are to guide 
the eye, and the dashed line is described in the text. 



value of the electron spin polarization V . On the other 
hand, for v^\ and low T, the well resonance W (Figs.|l|a 
and Rp) can also be much broader than this same fit 
(dotted lines). 

Figure shows the Ks(y) and the V^iy) of the W 
resonance near v — \^ for three different temperatures. 
As the temperature is lowered (Figs.0a-p|c), the sharp 
peak in Ksiy) evolves into a "tilted plateau". Figure 
0b also contains the Kgiy) data points reported pre- 
viously 0. While the new data are consistent with 
the earlier measurements, probing Kgiv) on this finer 
scale reveals a small region on both sides of v—\ where 
Ks{v)k.Ks{v=\)Iv (dashed line in Fig.|c). This tilted 
plateau is incompatible with the expression for Ks{y) 
derived previously l^], which had assumed delocalized 
quasiparticles |27[] . The existence of the tilted plateau 
is a natural consequence of the localization of the quasi- 
particles along the plane of the quantum well, such that 
the nuclei responsible for the W resonance see fully 
polarized electrons (7^=1), as if i^=l "locally", even 
though i/^l "globally". More precisely, for nuclei at 
'Ri={X',Y',Z'=0) in the center of the quantum well, 
the local KsmtCR-i) is directly proportional to the z- 
component of the local electron spin magnetization den- 
sity, MzCRi) p5| , which is in turn proportional to the 
product of the electron number density and the spin 
polarization, i.e., M^(Ri)oc |(?!)(Ri)P VCRi) |||. If the 
quasiholes (or quasiparticles) introduced into the system 
by going to v — 1 — e (or ly = 1 + e) are localized, then, in 
order to keep the total number of electrons fixed, |(/)(Ri) p 
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FIG. 3. KsiT) (filled symbols) and r^{T) (open sym- 
bols) for several filling factors 0.949<i^<1.045. Lines through 
Ks{T) are to guide the eye. Lines through r^o{T) are fits sug- 
gested by a simple model for Skyrmion dynamics (see text). 



must increase (or decrease) far from these charged exci- 
tations, which produces the observed tilt in the plateau 
near iy=l. This is just the same effect as the rise (or 
fall) of the water level in a pool induced by placing solid 
dipsticks (or hollow capillary tubes) into it to create lo- 
calized density minima (or maxima) , which correspond to 
quasiholes (or quasiparticles). In this analogy, changing 
the filling factor by adjusting i?totcos6' varies the number 
of dipsticks or capillary tubes in the pool and the Knight 
shift is given by the water level. Clearly, the fact that 
NMR is a local probe has become important. 

Taking different slices through the (v, T) plane pro- 
vides additional insights. Figure || shows Ks{T) and 
r^(T) for several filling factors near v—l. For lyj^l, 
lowering the temperature causes Tw{T) first to increase 
and then to drop, in stark contrast to the temperature 
independence of the well linewidth at i/=l. The non- 
monotonic temperature dependence (Fig.H) is consistent 
with the evolution of the W resonance from motionally- 
narrowed to frozen as the temperature is lowered. Quali- 
tatively similar trends were uncovered in earlier measure- 
ments at I' < 1/3 p3]. Both cases are rather unusual ex- 
amples of motional narrowing phenomena in NMR, since 
the nuclei are fixed in the lattice at such low T. Instead, 
the motion is that of delocalized spin-reversed quasipar- 
ticles, that results in fluctuations of the local hyperfine 
field JSf (Ri) at each nuclear site Ri. The shape of the 
resonance is sensitive to 6(Ri) = T(Ri)(5i?|(Ri)*7, where 
T(Ri) is the characteristic time scale of the fiuctuations, 
and '7 is the nuclear gyromagnetic ratio |25|. As T is 
lowered, the W resonance goes from the motionally nar- 
rowed limit (at high T, 8(Ri) << 1) to the "interme- 



diate limit" (at T near T^,,,, where T^{T^ax)=TT" ^ 
e(Ri) -- 1) and then to the "frozen hmit" (at low T, 
6(Ri) >> 1). Figures H and | show that all three hmits 
are experimentally accessible near v=\. 

Qualitatively, the |(5z^|-dependence of r^(T) (Fig. ||) 
appears to rule out several candidate mechanisms for 
the localization. For example, if individual Skyrmions 
were strongly pinned by a distribution of traps, then, as 
5v increases, quasiparticles should "start" to localize at 
the same temperature, but "finish" at lower and lower 
temperatures. In contrast, as 5v increases, the whole 
r^(r) peak shifts to lower T (Fig. |); e.g., at T w 1.5 
K, all Skyrmions can appear to be either localized (at 
\5v\ = 0.031) or delocalized (at \5v\ = 0.051). We con- 
clude that a collective process is required to explain the 
trends in Figs, g and y. This process, however, does 
not appear to be the "melting" of a classical Skyrmion 
crystal, since the classical melting temperature should 
increase as the crystal density (and the bond energy) 
increases. Two other collective mechanisms are qualita- 
tively consistent with the data. In the first scenario, the 
data are explained by the thermally assisted melting of a 
quantum Skyrmion crystal, which is approaching a quan- 
tum melting transition [g8| at some \5u\ > 0.05. In the 
second, the delocalization is due to the the"depinning" 
of a Skyrmion crystal, since the soft (stiff) bonds of the 
crystal at low (high) \5v\ may easily (not easily) stretch 
to match the disorder potential, resulting in a high (low) 
depinning temperature J29| . 

Additional, quantitative information about the local- 
ization process p6[ may be obtained by simulating the 
effect of Skyrmion dynamics on the OPNMR spectra 
@|6|. For example, to explain the r^(T, \5v\ = 0.031) 
data, we require X « 3 even as —^1 from below; the 
Skyrmion still has a "large" spin even as it takes ~ 20/iS 
to travel over the inter-Skyrmion spacing. Furthermore, 
in simulations that assume spatially-uniform local field 
fluctuations, we obtain F™"^ values that grow with in- 
creasing I (5 1' I, as seen in earlier measurements |2j]. In 
contrast. Fig. ^ shows that F™"^ is essentially constant 
for 0.031 < \6v\ < 0.051. Apparently, the actual localiza- 
tion process causes local field fluctuations that become 
more inhomogeneous as |(5i^| increases pi| . 

In the frozen limit, the fit previously used for the mo- 
tionally narrowed OPNMR spectrum |2^j2J] is replaced 
by the more general formula: 
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where the sum runs over all nuclei in the center of 
the quantum well, and Ksiat(R,i) oc Mz{R.i), which in 
turn reflects the shape of individual quasiparticles and 
their spatial arrangement along the plane of the quantum 
well. Since conventional Skyrmions are expected to form 
a square lattice, we consider this possibility first. 
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FIG. 4. (a) Expected radial dependence (in units of the 
magnetic length lo) of M'^{r), for the KO, K3, and K3* 
Skyrmions described in the text. For |5!/|=0.05, gray scale im- 
ages of Mf''"(Ri) (Black=-1; White=-|-1) are shown within 
the unit cell of a square lattice of either (b) KO, (c) K3, or 
(d) K3* Skyrmions. Also shown (in b-d) are black contour 
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Several theoretical approaches have been used to cal- 
culate M^(Ri) for a single Skyrmion excitation of the 
v = 1 ground state. Figure Ha shows the typical ra- 
dial dependence of the dimensionless M'^ (r) expected H] 
for Skyrmions with K=Q ("KO") and K=i ("K3") re- 
versed spins. The KO Skyrmion corresponds to the ordi- 
nary Laughlin quasiparticle, while both theory and ex- 
periment suggest that the K3 Skyrmion is energetically 
preferred for typical experimental conditions. Also shown 
in Fig. ^ is an ad hoc hybrid between the two ("K3*"), 
that has both the tail of KO and the three reversed spins 
of K3. 

We approximate MS£i{) for N localized Skyrmions 
with Mf*"(Ri) (Fig. |b-|d). Using Eqn. (|I|), we sim- 
ulate the frozen limit OPNMR spectrum at v=\ + 5v, 
for various Skyrmion shapes. Figure g] shows [Kg^™-{i'), 
^w^i^)] extracted from these simulations, which may 
be quantitatively compared to the low temperature data 
of Fig. ||c. We can rule out the standard model of a 
square lattice of conventional K3 Skyrmions over most of 
the plateau (see Fig. ||a and gb) ; using smaller conven- 
tional Skyrmions (K2, Kl) helps, but not enough. In- 
stead, the data are in much better agreement with simu- 
lations assuming a square lattice of either KO or ad-hoc 
K3* Skyrmions. Apparently, the existence of the tilted 
plateau requires that M^""(Ri) ~ 1 over a large frac- 
tion of the area between the quasiparticles, as in Figs. 
0b and Hd. This conclusion is not sensitive to the de- 
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FIG. 5. Open symbols are [Ks{i'), Tw{v)] data from Fig. 
2c. Filled points are [K§^'^{u), rS''"(i/)] extracted from sim- 
ulations described in the text. The points for each Skyrmion 
type are joined by lines. 



tails of our simulation over the range (|(5i^| < 0.05) of 
the observed plateau (e.g., changing to a triangular lat- 
tice, or including small disorder in Skyrmion locations). 
If the localized state is a 2D lattice of quasiparticles, it 
appears that they are either Laughlin quasiparticles (KO) 
or Skyrmions with very short tails, like the ad-hoc K3*. 

Alternatively, the localized state may involve "clumps" 
of Skyrmions, which result in large Skyrmion- free regions. 
This may happen, for example, if the disorder poten- 
tial favors large length scale density fluctuations, subject 
to the constraint of a collective localization process (see 
Figs. H and |^). More theoretical work is required before 
this picture can be compared to the data. 

For a 2D lattice, the Skyrmion shape consistent with 
the data is surprising, since the energetics of a single 
Skyrmion state would favor the conventional K3 over a 
shorter tailed Skyrmion (like K3*). However, this pref- 
erence may not be the same in a many Skyrmion state, 
where the energetics are more complicated. For exam- 
ple, the energy of a single Skyrmion is independent of the 
phase angle (p which defines the global orientation of the 
XY spin components. In the crystalline phase, however, 
interactions are generally expected to lead to preferred 
values for the relative phase angle {(j)i-<j)j) between the 
Skyrmions at sites (i,j) Q. Recent theoretical studies 
of a square lattice of K3 Skyrmions suggested that it was 
equivalent to the superconducting phase (SC) of the bo- 
son Hubbard model ||l^. We speculate that the shorter 
tail of the K3* Skyrmion would allow 0i,0j to remain un- 
correlated, thereby fixing the number of reversed spins on 
each site Ki to an integer, which corresponds to the Mott 
insulating phase (MI) of the boson Hubbard model [^ . 
If the total energy of this MI phase is lower than the SC 
phase, the system prefers the K3* Skyrmion. 
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